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Abstract 

The currently available set of hadron abundances at the SPS for central S 
+ Au(W,Pb) collisions is compared to predictions from a scenario assuming 
local thermal and hadrochemical equilibrium. The data are consistent with 
a freeze-out temperature T = 160 - 170 MeV. Spectra are consistent with 
this temperature range and a moderate transverse expansion. The freeze-out 
points at the AGS and SPS are found to be close to the phase boundary 
between a hadron gas and an ideal quark-gluon phase. 
Studies of ultra-relativistic nucleus-nucleus collisions in fixed target experiments at the 
BNL AGS and CERN SPS investigate hadronic matter at extreme density. Nuclei at the 
AGS are found to stop each other completely in the cm. frame, while at the SPS an onset 
of transparency is observed for S+A collisions [l[|. It has been predicted [@] that the longest 
lived systems with high density are produced somewhere between AGS and SPS energies. 
Hadronic cascade models yield densities in the interior of the colliding nuclei up to ten times 
normal nuclear matter density ||. An intriguing possibility is that the phase boundary 
to quark matter is crossed in these collisions. For zero net baryon density the transition 



temperature resulting from numerical simulations of QCD on a lattice is 150±10 MeV A 
first order transition can be constructed for finite baryon density or non-zero up/down quark 
chemical potential between an ideal hadron gas and an ideal quark-gluon gas and the result 
will be shown below. We will first address the question to what extent one is allowed to 
actually talk about a phase in the thermodynamic sense. The question is whether the time 
scale is long enough and/or the collisions are frequent enough for the system to equilibrate 
before it breaks up ( "freezes out" ) into the final state hadrons. Present data are now detailed 
enough to address whether the system is in equilibrium at freeze-out. 

Early data on global observables were seen || as indicative of a system with a temperature 
in the vicinity of 150 MeV, in local equilibrium but expanding longitudinally. Since then 
an extensive set of data on hadronic abundances and spectral distributions has become 
available for silicon and sulphur induced reactions and the hypothesis of thermal as well as 
hadrochemical equilibrium has been addressed by several authors |6Hl2|. In particular, in 
pTjj] an impressive and complete survey is given of theoretical techiques and their applicability 
to interpret data. For the AGS we have shown recently |TT| that the complete set of hadron 
abundances is consistent with a system in equilibrium at a temperature in the range 120 
- 140 MeV, a baryon chemical potential of 540 MeV, and strangeness in equilibrium with 
up/down flavors. The system appears to be in local equilibrium with an overall longitudinal 
and transverse expansion with average velocities of (/?;) = 0.52 and (/3 t ) = 0.39 - 0.33. Here 
we will use the same approach to test whether data at SPS energies are consistent with 
hadrochemical equilibrium as well. 

We use for the present analysis the complete set of data now available from the differ- 
ent experiments. Previous studies often concentrated largely on the abundances of strange 
and multiply strange baryons, in part by choice and in part because some of the relevant 
nonstrange hadron and meson abundances were not available at the time. This approach 
yields significant discrepancies between model predictions and data when the most abun- 
dantly produced particles, pions, are considered, as was already noted by Redlich at al. 
0. Further, we will use data integrated over the maximum available range of rapidity and 
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transverse momentum since flow effects can severely distort relative hadron abundances at a 
given rapidity and transverse momentum. Similar to |||7| we will investigate only the case of 
complete strangeness saturation instead of allowing an additional free parameter to govern 



abundances of strange hadrons as in p|- [T0|JT^] . In addition to comparing relative hadron 
abundances to the data we will also compare the absolute pion density from the thermal 
model to an estimate obtained from the experimental data. 

The basic assumption of our thermal model is that in every local restframe the system is 
described by a grand canonical ensemble of fermions and bosons in equilibrium at (freeze-out) 
temperature T. For an infinite volume the particle number densities are given as integrals 
over particle momentum p: 

Pl 2tt 2 J exp[(E l - f i b B t -fi s S i )/T]±l U 
where gi is the spin-isospin degeneracy of particle i, E iy E>i and Si are its total energy 
in the local rest frame, baryon number and strangeness, and fib and fi s are the baryon 
and strangeness chemical potentials (unless otherwise noted h =c=l). Energy density e°, 
pressure P®, and entropy density of for a given species are obtained by evaluating equation [I] 
above with an additional factor of Ei, p 2 /(3Ei), and (p 2 /(3Ei) — + Ej)/T in the integrand, 
respectively. 

For a system of finite size the integrand in equation (1) has to be multiplied by a correc- 
tion factor ||13|| . For an estimate of this correction we assume a spherical volume with radius 
R giving a correction factor 

37T 1 

f = 1 ~ 4p~R + J^RY' (2) 
To approximately account for the volume taken up by baryons we apply an excluded volume 
correction to the partition function, 

In 7° 

\ nZi = ^ jr (3) 

where Vj is the volume occupied by an individual baryon and the sum extends over all 
baryons; we use a sharp sphere volume with radius 0.8 fm for all baryons. The excluded vol- 



ume correction for the pressure takes the same form as eq. |3|; to correct the entropy, particle, 
and energy densities the appropriate derivatives of the correction factor with respect to tem- 
perature and/or chemical potential have to be taken into account. This simple correction is 
appropriate for J2j VjP°j < 0.5, valid for the area of interest in this paper. More general (and 



much more involved) procedures are discussed in [H|. The finite size and excluded volume 
corrections are sizeable and affect the absolute densities and pressures, but ratios of particle 
yields or quantities such as the entropy/baryon are affected to a lesser extent. 

For the comparison to SPS heavy ion data we include in the thermal model all known 
fl5fl baryons and mesons up to a mass of 2 and 1.5 GeV, respectively. We have checked 
that for T < 180 MeV higher mass mesons and baryons do not play an important role and 
lead to corrections of less than 5 % in the particle densities. The correction for feeding and 



decay is performed as in [T]J using all known branching ratios |T5[ and symmetry and phase 
space arguments for unknown branching ratios. Unless an explicit number is quoted by the 
experiment, we assume a 50 % particle identification efficiency following weak decays. 

The starting point for our thermal model calculations is to determine, for each tem- 
perature, that baryon chemical potential with which are best described experimental data 



reflecting the proton to pion ratio (first four rows in Table |). As in |TT[| the strangeness 
chemical potential is fixed by the strangeness neutrality condition. Best overall agreement 
with all data is obtained with T = 160 - 170 MeV, as shown in Table | where all currently 
available experimental data on particle ratios measured in central S + Au(W,Pb) collisions 
are compared to predictions of the thermal model. A graphic illustration of what drives 
the temperature in our freeze-out analysis can be seen in Figure [l|, where predictions for 
those particle ratios that show, after our choice of fib, the largest temperature sensitivity are 
plotted vs T. Typical changes are about a factor 30 over the temperature range considered 
and good overall agreement with the data for these temperature sensitive ratios is obtained 
for T = 160 - 170 MeV and corresponding baryon chemical potentials of 170 - 180 MeV, 
respectively. Since we assume strangeness equilibration all other particle ratios, including 
those for strange and multi-strange hadrons, are then fixed. Inspection of Table | shows that 



with these two parameters one can obtain a surprisingly good description of the available 
experimental data. 

Most measured particle ratios are reproduced well within the errors, especially keeping 
in mind systematic errors and systematic acceptance effects as displayed e.g. in Table | by 
data for the same ratio from different experiments. Of the 27 measured ratios spanning a 
range of a factor 400 just three differ beyond statistical errors by 40-50 % and the only more 
serious discrepancy is the ratio (Q + + Q~)/(E + + E~) which appears to be about a factor 2 
different. Note, however, that there is no systematic indication that a separate parameter is 
needed to control strangeness abundance. At the same time, there are significant differences 
(up to a factor of four) in the ratios for the heavy ion reaction considered here as compared 
to nucleon nucleon data |3(| . 



The freeze-out temperature resulting from our analysis is considerably lower than that 
proposed by || and |1(J to describe data. Although in the latter a range of temperatures 
is considered, the finally proposed chemical freeze-out temperature is 190 MeV causing, 
as the authors note, a significant underprediction in the pion abundance as compared to 
data. The lower temperature is reassuring, because temperatures around 200 MeV will lead 
to absolute pion densities of 0.6-0.7/fm 3 after excluded volume correction, i.e. even pions 
start to overlap significantly (as noted also in flOfl). Furthermore, such a high freeze-out 



temperature would imply that freeze-out takes place well in the quark-gluon plasma region 
of the phase diagram (see below). In fact, experimental data on two-pion interferometry can 
be used to obtain an estimate of the pion density at freeze-out. We begin by noting that 
experimental tt + tt + and tt~tt~ correlations are rather well described by a simulation using 
the cascade code RQMD |31|. Next we inspect the space-time history of particles in RQMD. 
Pions freeze out over an extended time in RQMD during which the source is expanding. To 
give a typical volume we quote here the size of the system at the average pion freeze-out time 
t = 14 fm/c. At mid-rapidity the distribution of positions transverse to the beam is Gaussian 
with standard deviation a x = a y = 4.0 fm, the longitudinal distribution for low transverse 
momenta (pt < 0.1 GeV/c) has a standard deviation of a z = 7.7 fm. This corresponds to a 



volume V = (27r) 3 ' 2 a x a y a z = 1940 fm 3 . The rapidity integrated multiplicity for negatively 
charged particles for central S+Au collisions has recently been reported by NA35 |l7j and 
from their data we estimate that the total pion multiplicity in these collisions is about 
510. Using the volume estimate just given one obtains a typical pion density at freeze-out 
p n ~ 0.27/fm 3 . This is rather close to the value of 0.30/fm 3 obtained in the thermal model 
for T = 160 MeV and the above parameters. Other volume estimates, e.g. obtained by 
using the time averaged variances or the variances of the time integrated distributions, are 
larger, leading to lower freeze-out density estimates. 

The entropy at freeze-out, when evaluated in the thermal model, is found to be large. 
For T = 160 and 170 MeV we obtain for the entropy per net baryon, which is equivalent to 
the entropy produced per incident nucleon, values of S/(B-B) = 45.4 and 36.7, respectively. 

The strange chemical potentials obtained for T = 160 (170) MeV and //& = 170 (180) 
MeV are /i s = 38.0 (47.0) MeV. On the quark level this implies a strange quark chemical 
potential of fj, qs = — /i s = 18.6 (13.0) MeV, much smaller than the strange quark mass. 
Obviously, for every temperature T there is one value of //& which, together with strangeness 
conservation, will yield fj, qs = 0. In the present case the pion to nucleon ratio happens to be 
such that we are close to this situation. This situation arises in a purely hadronic picture and 
is not related to whether or not the system is in the quark- gluon plasma phase at freeze-out. 

In order to test whether the transverse momentum spectra of the various hadrons are 
consistent with a temperature of 160 MeV and one common transverse flow velocity we 
use, as in our analysis for the AGS JTTJ, the formalism developed in [|32[| . To minimize 
systematic uncertainties, we choose spectra of different particle species measured in the 
same experiment and close to mid-rapidity and we have made the comparison for three sets: 



pion, kaon, proton and deuteron spectra from NA44 [33|, where a similar fit is also shown; 



spectra of kaons, lambdas and cascades from WA85 p8|j24| shown together with the fit in 
Fig. U the f] to pion ratio from WA80 [ 2^| . All spectra are consistent with a temperature of 
160 MeV and, for a linear dependence of the flow velocity on the radius of the system (a = 
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1), an average transverse expansion velocity of {(3 t ) = 0.27. Consistent results were reported 
by p2[ for spectra from S+S collisions. This value is somewhat smaller than the range of 



0.33 - 0.39 found at AGS energies 11 



To put the above determined freeze-out temperature and baryon chemical potential into 
perspective we have calculated the phase boundary between the hadron resonance gas and 
the quark-gluon plasma by equating the chemical potentials and the pressure in the hadronic 
phase, with those of an idealized phase of massless u,d quarks, s quarks with mass of 150 
MeV, gluons, and a bag constant of B = 262 MeV/fm 3 . The resulting phase diagram is 
shown in Figure ^| along with the latent heat and baryon densities at the transition line. 
Interactions among the hadrons, which are neglected in our approach, are not expected to 
change this phase boundary by much except possibly at high baryon density (fib > 1 



GeV), i.e. far away from the freeze-out region. The freeze-out points determined from the 



present analysis and from that for AGS data |TT| are shown by the filled circles in this 
diagram. Note that, with this scenario, the system at freeze-out, i.e. after expansion and 
cooling, is close to the phase boundary at both AGS and CERN energies. 
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TABLES 

TABLE I. Particle ratios calculated in a thermal model for temperatures of 160 and 170 MeV, 
baryon chemical potential /ib of 170 and 180 MeV and strangeness chemical potential fi s of 38.0 and 
47.0 MeV, in comparison to experimental data (with statistical errors in parentheses) for central 
collisions of 200 A GeV/c S + Au(W,Pb). For experimental data a p t range is quoted when the 
lower limit is significantly larger than zero. 
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FIGURES 

FIG. 1. Ratios of baryon abundances that show a strong temperature dependence for a fixed 
pion to nucleon ratio as indicated in Table |. See text for details. 

FIG. 2. Kaon, Lambda and Cascade spectra from WA85 [28,24] compared to thermal model 
calculations with T = 160 MeV and an average transverse expansion velocity of {(3t) = 0.27 (solid 
lines). 

FIG. 3. Phase boundary between a hadron gas and a quark-gluon plasma (top) as function 
of temperature and baryon chemical potential together with the freeze-out points for Si(S) + 
Au(W,Pb) collisions at AGS [11] and SPS (present paper) energies. Latent heat of the phase 
transition (middle) and baryon density in the hadron and quark-gluon phase at the phase boundary 
(bottom). 
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